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Summary. Calculating an exact self energy for ab initio transport calculations 
relevant to Molecular Electronics can be troublesome. Errors or insufficient approxi- 
mations made at this step are often the reason why many molecular transport studies 
become inconclusive. We propose a simple and efficient approximation scheme, that 
follows from interpreting the self energy as an absorbing boundary condition of an 
effective Schrodinger equation. In order to explain the basic idea, a broad intro- 
duction into the physics incorporated in these self energies is given. The method is 
further illustrated using a tight binding wire as a toy model. Finally, also more re- 
alistic applications for transport calculations based on the density functional theory 
are included. 



1 Introduction 

The most important driving force in the research field of Molecular Electronics 
are prospects on technological applications - whence the name - in entirely 
new realms of system parameters [1]. The development of these new tech- 
nologies also requires serious progress in several disciplines of fundamental 
sciences including both, theory and experiment. One of the major theoretical 
challenges is the quantitative description of transport through a molecule with 
a given contact geometry [2, 3, 4, 5, 6]. 

In order to appreciate the caliber of the problem, recall that describing 
transport requires to keep track of two aspects of physical reality, each by 
itself posing a task of considerable difficulty. Needed are a) a good know- 
ledge of molecular states, i. e. energy levels and orbitals, which is not easy to 
obtain, since they experience a strong influence by Coulomb interactions on 
the molecule, and b) a thorough understanding of the hybridization of these 
orbitals with the electronic lead states, so as to predict the broadening, i. e. 
the "life time" , of molecular energy levels. This seriously complicates matters 
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for ab initio calculations, because inevitably a macroscopic number of degrees 
of freedom is involved. We arc facing here a classical dilemma: a single one of 
the two problems - interactions on the molecule and the macroscopic number 
of lead atoms - by itself can be dealt with reasonably well, but only at the 
expense of applying methods that exclude a simultaneous solution of the other 
problem. In a sense, we find ourselves in a situation not unlike Ulysses, when 
he was trying to pass by Scylla and Charybdis [7] . 

In this paper we present a method, that simplifies b), i. e. including macro- 
scopic electrodes into ab initio calculations. The incarnation, that we put for- 
ward in this communication, operates in those instances where the calculation 
of transport coefficients builds upon a formalism in terms of Green's functions. 
The basic idea developed here, however, is much more general than that and 
may also be of use for example in transport calculations based on the density 
matrix renormalization group [8] . 

A typical example of a Green's function based transport theory met 
in cases, where the quasi-particles are effectively non-interacting, is the 
Landauer-Biittiker approach to transport [9, 10]. In it the conductance (in 
units e^/h) is expressed as a transmission at the Fermi energy, g=T{E-p). 
Explicitly, T{E) has a representation 

T{E) =tr GFiG'^r, (1) 

which may be derived using elementary scattering theory [11, 12, 3], the 
Keldysh technique [13] or the Kubo formula [14, 15], in principle. 

The "dressed" Green's func;tions, G, required in any of these approaches 
describe the propagation of particles with energy E on the molecule in the 
presence of the electrodes. The external leads, left and right, enter these func- 
tions by self energy contributions, ^f,,, one for every electrode. They relate 
G to the Green's function of the isolated molecule, G^ by 

G-\E) = G-\E)-S^{E) (2) 

and include all the effects of coupling to the left and right leads, Sm=Si:+S,. 
Also, they determine the level broadening F, ^r=i{^i .r^^] r) appearing in 
(2). ( Equation (2) should be understood as a family of matrix equations with 
resolvent operators G, G^i , parameterized by energy, E; Green's functions are 
actually the matrix elements of these operators, x, x')=(x|G(£')|x') and 
G«(^;,x,x')=(x|G«(£;)|x').) 

The calculation of the exact couplings, Si^r, usually is fairly troublesome. 
In the simplest case, when the electron interaction can be appropriately dealt 
with by an effective single particle model, the couplings take a structure 

S,=t,Gs-J„ K = C,r (3) 

where denote the two hopping matrices that connect the molecular junc- 
tion with the left and right electrodes [13]. The "surface" Green's function. 
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Gj-Cn debuting here describes the propagation of quasi-particles on the elec- 
trodes in the presence of the contact surface. Even in this situation calculating 

is not really easy. Complications arise since a) should include macro- 
scopic leads plus contact geometry and b) the hopping matrix t is (normally) 
not just of the nearest neighbor type, so the contact surface also involves 
sub-surface layers of electrode atoms, in general. 

The procedure to be proposed in this communication simplifies the con- 
ductance calculation by essentially eliminating the step of evaluating S.^ . It 
works after the molecule has been redefined. The "extended" molecule, esW , 
not only includes the original molecule but also pieces of the left and right 
electrodes: 

G-\E)^G-^-S,^. (4) 

Our key observation is the following: while the molecular conductance is cru- 
cially dependent on microscopic details incorporated in , it is completely 
indifferent towards details in S^^^ , ij eM includes sufficiently many electrode 
atoms. As a consequence, there is no need to use the exact self energy in 
order to obtain (in principle) exact results. ^ One can replace it by a simplistic 
model coupling of the type 

(x| S,^ {E) |x') ^ ir?(x) 5(x - x'), (5) 

where we have introduced a "local leakage function" r]{x). It is crucial to our 
method, that fine tuning ry(x) is obsolete once certain criteria to be specified 
in Sec. 2.3 are satisfied. 

The outline of this paper is as follows. In Sec. 2 we recapitulate in broad 
terms the physical effects, that are encoded in the self energy formalism. 
The concept of the extended molecule will emerge quite naturally from these 
considerations. They will also illuminate under what conditions (5) can be 
justified. In the following Sees. 3 and 4 we will present a series of model 
problems. In order to demonstrate the principle, we begin in Sec. 3 with 
a tight binding chain, for which numerical results can be compared against 
analytical solutions. To illustrate the usefulness in practically relevant cases, 
the conductance of di-thiophenyl is investigated in Sec. 4 using an approach to 
transport based on the density functional theory and the quantum chemistry 
package TURBOMOLE [16, 17, 18]. This study is also intended to reveal the 
limits of the ansatz (5). 



^ The words "exact result" have in the present context a restricted meaning: they 
refer to the exact solution of the single particle scattering problem, that can 
be stated once the Kohn-Sham orbitals and energies are given. Under which 
conditions - if at all - scattering theory based on (ground state) Kohn-Sham 
orbitals could give an exact description of the full many body problem, this is 
an important question which, however, goes well beyond the scope of the present 
article. 
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2 Basic Physics of Self Energies 

In this section we will give a more precise definition and a justification of 
the procedure (5) for constructing a self energy, which is based on physical 
arguments. In order to explain the logic, wc first recall on a qualitative level, 
which physical information is carried by the original self energy, ■ We shall 
illustrate then, how this information is transfered into Ggsw by reformulating 
the problem to calculate G in terms of an extended molecule. If enough metal 
atoms have been included in eM , the "information transfer" will be complete. 
Then the remaining information in the self energy of eM is trivial, i. e. it 
is no longer molecule specific. Therefore, SgM is apt to simple approximations 
like (5). 

2.1 Self Energy of the Molecule Sm 

The self energy, , that appears in (2), 

S^{E)=G-\E)-G-\E) 

has two qualitatively different effects which are incorporated into its hermitian 
and anti-hermitian constituents. 

Hermitian Constituent 

The eigenvalues of GeM for the isolated molecule are real numbers. Due to the 
hermitian piece, (Jiffif =(i75if +Z']^)/2, those eigenvalues undergo a shift, Ae^, 
when the molecule is coupled to the electrodes. In the case of weak electron- 
electron interaction this simple "renormalization" of excitation energies is all 
that can happen. However, if the interaction is strong, such that the electrons 
are highly correlated, additional and qualitatively different effects can occur. 
A most prominent representative is the Kondo effect, observation of which has 
been reported in various recent experiments [19, 20, 21]. It manifests itself in 
the spectral function of the coupled molecule 

A{E) = (V27r) Tr{G{E) - G^ (£)), (6) 

which measures the number of molecular excitations with a certain energy, 
roughly speaking [22]. Kondo-physics is signalized by an additional peak in 
A{E), the "Abrikosov-Suhl" -resonance, which sits right at the Fermi-energy of 
the leads [23]. This resonance is a collective phenomenon involving electrons 
from the leads and the molecule; it cannot be understood as a renormalization 
of a molecular energy level alone. 

Even in the absence of strong correlation effects, the shift of molecular 
excitation levels, Ae^, can have very important consequences for the inter- 
pretation of experimental findings. The presence of the metal electrodes can 
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help screening the interaction of electrons on the molecule. As a consequence, 
the energy difference between the lowest unoccupied molecular energy level 
(LUMO) and the highest occupied level (HOMO) will generally shrink. In 
the extreme case, where the LUMO falls below (or the HOMO above) the 
Fermi energy of the electrodes, charge will flow onto the molecule such that 
the molecular junction becomes partially polarized even at equilibrium con- 
ditions. 

Anti-Hermitian Constituent: Exponential Time Evolution and 

Reservoirs 

The hcrmitian piece of Z'af is basically "just" a modification of the (effec- 
tive) Hamiltonian. By contrast, the anti-hermitian piece of the self energy, 
{S^ — sl^)/2, introduces a qualitatively new aspect, because it gives rise to 
an imaginary component, ij^, of the eigenvalues of G. It gives the molecular 
levels, u, a finite lifetime reflecting a simple physical fact: an initial excitation, 
localized at time t=0 on the molecule, can fade away to be absorbed by the 
leads, ultimately. 

Let us discuss how excitations pass away in more detail so as to see, why 
the electrodes and the thermodynamic limit are important ingredients in un- 
derstanding the self energy. We begin by noting, that Green's functions can de- 
scribe a time evolution of the physical system. Therefore, the relaxation rates, 
7i^, also have a straightforward interpretation in time space. Assume, that the 
molecular junction is prepared in an initial state such that the molecule has 
an excess charge. Then, the rates 7^ describe an exponential decay in time, 
exp(— 7i,t), exhibited by each contribution to this charge made from a certain 
molecular level, i/. 

Now, the exponential dependence exposed here is implied to be valid at 
all times including in particular the asymptotic regime t — > 00. This means, 
that the charge is really swallowed by the electrodes, it never returns to the 
molecule and only for this reason the relaxation process can ever become 
complete. In other words, the electrodes act like thermodynamic baths or 
reservoirs. They destroy information about the initial state in the sense that 
the return time of a signal, i. c. electrons, from the reservoirs is infinite. 

As usual, a truly diverging time scale can be realized only with infinitely 
many degrees of freedom; otherwise return paths (c. g. of electrons) exist with 
an overall weight that is not vanishing. In this infinite dimensional Hilbcrt 
space the time evolution is unitary, of course. The (anti-hermitian part of 
the) self energy pops up as a consequence of projecting the full time evolution 
down to the subspace of the molecule, W , which then can no longer remain 
unitary. The principle encountered here is well known in the general theory 
of non-equilibrium phenomena [24]. 
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Self Energy and Transport 

Clearly, the decay rates 7y must be closely related to transport properties, 
because they govern the time evolution of charge exchange between molecules 
and leads. Note, however, that the self energy of the Green's function contains 
information only about the total loss rate, 



due to leakage. It does not necessarily keep track of the rates Z'f ^ separately, 
that describe the exchange with the individual leads, left or right. This latter 
piece of information is important for the transport characteristics, as can be 
seen e.g. in the Landauer-Biittikcr formula (1). In general, it cannot be recon- 
structed from the G{E, x, x') alone, without making further assumptions (e.g. 
that {x\Em |x') is block-diagonal with the two diagonal entries resemblingi^r^f , 
separately) . 

In order to illustrate the significance of the level shifts Ae„ and level broad- 
enings for the transport problem, we consider now a situation typical of 
experiments on molecular conductance. We focus on the case of weakly inter- 
acting electrons and call 6hi the energy gap between the HOMO of the isolated 
molecule, Ch and its LUMO, ef. Sh\=ei—eh- In typical transport experiments 
one has a situation, where 5hi > lel^. At the same time, the experimentally 
measured values of the conductance, g, of the molecule only very rarely exceed 
0.1. Both observations taken together give a strong indication that for this 
type of experiments the level broadening of HOMO and LUMO, 7h.i, is well 
below the level separation, 7h,i <C ^hi- Roughly speaking, the conductance (1) 
will then be given by a superposition of two Lorenzians, 



with a Fermi energy of the metal, E-p, situated in between the values of HOMO 
and LUMO after coupling, eh.l = eh,i + ^eh,i- 

We add a remark regarding uncertainties in theoretical predictions of level 
positions and their broadenings. Inaccuracies in calculating absolute values of 
the level positions tend to induce a shift of the transmission curve, but do not 
normally change their structure - unless molecular levels happen to cross the 
Fermi energy of the electrodes, of course. Often, the shift is very similar for 
all energy levels involved, and therefore it can be partially eliminated when 
the transmission is plotted over E — Ep . 

Inaccuracies in the level broadening are more severe, since their error turns 
out to be of the order of unity. The value of the conductivity off resonance is 
determined by 7^ j, and so a quantitative calculation of g under these condi- 
tions is very difficult. The source of this error and the question how it can be 
overcome became a very active field of research, recently [2, 3, 4, 5, 6]. 
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Relation to the Renormalization Group Method 

In this section, we describe the physics incorporated in from the point of 
view of an hypothetical renormalization group method. This is to say, that we 
investigate how evolves when we build up the molecular junction gradually 
step by step, attaching more and more electrode atoms. The idea is in a 
spirit similar to the density matrix renormalization group [25] . The flow thus 
induced will be smooth unless the molecule becomes strongly distorted, which 
could signalize for example dissociation or ionization. 

In order to illustrate this evolutionary process, we have performed cal- 
culations based on the density functional theory (DFT) using the standard 
functional BP86 [26, 27]. DFT provides us with an effective single particle 
Hamiltonian, Hj^,^, with eigenvalues and corresponding eigenstates Our 
interest is in how the eigenvalues and eigenfunctions change when we include 
a gradually increasing number of electrode atoms, , in our model system. 

The result of this procedure has been depicted in Fig. 1 for the case of the 
molecule di-thiophenyl (See Fig. 4 for the detailed atomic structure.) Every 
eigenfunction is represented by a data point (ejy, A^). The integrated ampli- 
tude is defined as 

Ai, = Tr« \v){v\ 

where the Tr,^ is over the projected segment of the Hilbert space that is 
associated with the molecular degrees of freedom. Our calculation is performed 
using a local basis set \Xi) (TZVPP [28]), with basis functions labeled by 
atomic positions, X, and orbital quantum numbers, £. When evaluated in this 
basis, the Tisi^ is a sum over all basis states that belong to molecular (i. e. 
non-Au) atoms. 

Fig. 1 shows, how the six orbitals of the molecule in the gas-phase shift 
and hybridize with TVj; increasing from 0, 1,5, 14,30,71. For illustration, we 
have also given the Kohn-Sham wavefunctions of two representative states in 
Fig. 2. 

The overall plot clearly shows, that the original orbitals survive the cou- 
pling to the electrodes and therefore contribute as resonances to the transport 
characteristics. The initial evolution at small N.^ is not very smooth, which is 
because a) attaching the first few Au-atoms cannot be considered a very small 
perturbation to the molecular system and b) at "magical" numbers of atoms, 
e.g. Njr=5, the electrode configuration is particularly stable. These "stabil- 
ity" islands are interesting in themselves but for our present purpose they 
deliver parasitic side effects, since they make it more difficult to extrapolate 
the overall flow. When keeping away from exceptional numbers, e. g. taking 
A^2=14, 30, 71, the evolution shows the expected smooth behavior. 

We have already mentioned, that the smooth (n'olution of single particle 
levels can also be perturbed, if a prominent molecular level happens to cross 
the Fermi energy. This can happen in a situation where the HOMO is relatively 
close to Ep. Then small fluctuations of the charge distribution, that occur 
due to the gradual appearence of "evanescent" modes, i.e. invading electrode 




Fig. 1. Di-thiophenyl between Au-pyramids (symmetric coupling, see Fig. 4). Left 
panel: flow of energy of different orbitals with increasing number, A'l , of electrode 
Au-atoms as indicated at the right hand side Eixis. (Orbitals at a),b),c),d) are shown 
in Fig. 2.) Each orbital is characterized by its energy and weight A„ on the molecule. 
The vertical bars near — 5eV mark the (center of the) HOMO-LUMO gap. The 
evolution of the six orbitals of the isolated molecule has been indicated by "world" 
lines. Right panel: evolution of charge, Q, accumulated on the molecule (including 
sulphor atoms). 



states with energies in between the prominent molecular orbitals, can lift the 
(designated) HOMO above Ep at certain electrode configurations. 

This is what is happening with Au55 electrodes, as can be seen from Fig. 3. 
In this case a molecular state (that turns out to be localized predominantly on 
the S-atoms) peaks above Ep and therefore is evacuated, leading to a decrease 
of the charge accumulated on the molecule by 1.2e. The fate of this prominent 
mode that has been expelled from the region of occupied energy levels is a fast 
decay with further increasing N.^ , because of its relatively strong coupling to 
the "invaders" , see Fig. 3. 

So far we have witnessed the transformation (or in some cases the decay) 
of the states, that occurs as a consequence of the hybridization of electrode 
and molecular orbitals. In a sense, this is the analogue of Fermi-liquid theory. 
The Kondo-effect, which in principle could appear for molecular systems that 
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Fig. 2. Plots of two orbitals, ij), for di-thiophenyl without electrodes (left) and 
with Ni,=<iQ Au-atoms attached (right). Orbital a) hybridizes only very little with 
the electrodes and is almost unaffected after coupling, b). By contrast, orbital c) 
becomes completely delocalized in the metal, see d). Grey (black) colors indicate 
regions with ip positive (negative). 

carry a spin, cannot be understood within this picture. This is because the 
Abrikosov-Suhl resonance is not a shifted molecular level. Instead, it is a 
collective phenomenon and generated by a large number of electrode states. 
Their energies reside inside HOMO-LUMO gap of the "dressed" molecule 
close to the Fermi-energy. This effect can be seen with DFT in principle, but 
the current approximations for the exchange correlation functional are too 
crude in order to capture it. A study with the exact, at present unknown 
functional should show satellites at Ey produced by lead states that merge 
with one another if the system size becomes large. If sufficiently many of them 
superimpose, a sharp peak, the Abrikosov-Suhl resonance, grows right within 
the HOMO-LUMO gap. A further typical characteristics of this collective 
effect is, that the resonance position is always the Fermi energy, irrespective 
of shifts in the molecular orbital energies, that might be induced by a gate, 
for instance. 

2.2 Extended Molecule and 

The self energy of the original molecule, , can contain a wealth of nontrivial 
information, it is not a quantity easy to calculate. This was the message of 
the preceeding section. However, the situation greatly simplifies, after the 
molecule has been redefined. Let us consider an extended molecule, eM , that 
comprises in addition to the original molecule also a "contact region" , i. e. a 
number N,^ of electrode atoms. The Green's function for an extended system, 
GeM , is related to the full Green's function G via a new self energy 



G-\E)^G-2{E)-S,,,{E). 



(8) 
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Fig. 3. Plot similar to previous Fig. 1 with A'^r=55 included. A right shift of the 
prominent orbitals is caused in the center panel because of additional, "evanescent" 
modes appearing. 



We give two reasons, why it is that S^^i is much easier to handle than . 

Imagine the extreme limit, in which far more electrons are sitting on the 
metal than on the molecule. Then the HOMO of the big system, en, is given 
by the Fermi energy of the metal, Ep, up to a small uncertainty, which is 
of the order of the HOMO-LUMO gap of eiW, Jhl- In a metal, this gap is 
inversely proportional to the number of metal electrons in the calculation, 
^HL oc 1 /A^eM , so that the uncertainty of the position of the Fermi energy 
with respect to eH,L can be made arbitrarily small. This is a very obvious 
advantage. 



Fig. 4. Atomic configurations of molecular junctions attached to pieces of the elec- 
trodes. Electrodes are modeled by pyramids, A'^i=14 Au-atoms each. Dark atoms are 
surface atoms used in self energy (29). Left: 2 Au-atom wire, Ns — 13 (Sec. 4.1). Cen- 
ter: di-thiophenyl with (stretched) S — Aui coupling; 2 sulphor, six carbon atoms, 
A'^=13 (Sec. 4.2). Right: di-thiophenyl with relaxed S — AU3 coupling, A'^^ll (Sec. 
4.2). 
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More importantly, the flow of the typical level broadening, 7h,l (related to 
Se!u)i that is driven by increasing the number of electrode atoms, N^, in the 
contact region, will lead us into a very tractable regime as we shall see now. 
It is only a fraction of the N.^ electrode atoms that is actually connected 
to the rest of the leads, the "outside world". Assume, that A^^ grows in such 
a fashion, that the number of these "surface atoms" , A^^ , does not change, so 
we build a quasi-one-dimensional wire. Then, increasing implies, that the 
fraction of the wavcfunction amplitude of extended orbitals sitting near the 
surface decays like N^/N^, so that 7h,l scales like 6ul- In good metals the 
ratio of both energies is of the order of the metallic conductance 

7h,l/&l ~ g > 1. (9) 

So, the level broadening of HOMO and LUMO of the extended molecule 
always exceeds their separation if the electrodes are made from a good metal. 
This is a situation, exactly opposite to the problematic one, 7h.i/^hi ^ 1, 
that we have encountered before in the context of (7). Summarizing, for the 
extended molecule, eM , the following hierarchy of inequalities holds true (10). 

Skl < 7h,l < 7h,i < Sh\- (10) 

The separation of energy scales implied by (10) is the prerequisite for the 
real gain that one makes when one turns to the extended molecule. The point 
is that the fine structure in the spectrum of GeM is on the scale of (5hl- The 
anti-hermitian constituent of SeM ^ 7h,Lj provides the smearing of this fine 
structure necessary in order to obtain smooth curves, e.g. for spectral and 
transmission functions. The details of this smearing have very little impact 
on the resulting curves, because the interesting structures live on energy scales 
Sh,i and 7h,i, which exceed ^hl and 7h,l by a parametrically large factor. 

2.3 SeM and Absorbing Boundary Conditions 

There is yet another, perhaps particularly intuitive way to understand the 
principal difference between 17,^ and SgM ■ It will serve as a motivation for 
the proposed approximation (5). 

Let us assume, we opted for an investigation of transport properties in the 
time domain, e. g. by propagating wavepackets. Then, we would study the 
time evolution of a wavepacket, that is localized at t=0 at some initial position 
on the; molecule. In particular, we can investigate how wavepackets leak out 
of the molecule into the contacts such that they gradually disappear. When 
performing such an investigation systematically in the presense of leads, one 
can in principle collect enough information in order to reconstruct the Fourier 
transform of the retarded Green's function, G{t). 

There is a condition on the observation time, T. In order to have an energy 
resolution 7;_,. we need T > 7; This comes for us with trouble, if the contact 
size maintained in our calculation is not sufficiently large. After some time. 
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the leakage will hit the outside walls of the contact, reflect back off them and 
finally, after the "dwell" time, tjj, it will rcarrivc at the molecule; we calculate 
Gf,M{t) instead of G{t). The energy resolution, that we achieve with such a 
calculation, can not exceed /i/td- The best that we can hope for is td ~ 
but only in a case, where the contact acts as a fully chaotic cavity. At longer 
times the signal from the decaying wavepacket will be superimposed by the 
cavity modes that describe how the wavepacket sloshes back and forth inside 
the electrodes. We shall demonstrate this effect in Sec. 3 looking at an explicit 
example. 

The salient point we wish to make is, that the minimum dwell time in the 

cavity should be so long that the wavepacket has enough time to evacuate 
before the molecule is being refilled again from the backscattered modes: 

TDlCr » 1. (11) 

There is a very elegant and powerful procedure that eliminates spurious cavity 
modes so that the condition (11) is always satisfied: one introduces absorbing 
boundary conditions (abc) in some regions of the cavity. These "surface" re- 
gions should swallow incoming signals, i.e. wave packet amplitude, and thus 
properly mimic escape to infinity. If the wavepackets, bouncing hence and 
forth inside the cavity, are completely eliminated before the return time td has 
passed, no trace of the finite cavity/electrode size will be left and G(t, x, x') 
can be reconstructed. This is exactly how adding the self energy S^^ works 
and nothing more than this is implied, if the contacts are sufficiently large. 
Therefore, we are entirely free to replace the exact boundary conditions E^m 
by any other ones, provided they absorb sufficiently fast (and do not disturb 
the immediate vicinity of the molecule-electrode junction). 

This is the idea underlying the step proposed in (5) . From what has been 
said above, it should have become clear that this ansatz is actually not just 
a good approximation, but it will give exact results, if the number of metal 
atoms is sufficiently large and the damping function has been chosen well 
enough. The examples given in Sec. 4 suggest, that a relatively small number 
of contact atoms ~ 10 — 20 can already give reasonable results. 

3 Toy Models 

In this section we are going to analyze two toy problems as test cases, namely 
the conductance of an L-sitc tight binding wire, clean and in the presence of 
an obstacle, that mimic;s a molecule. 

We will begin with a single channel wire and show, that the technique 
introduced with (7) delivers the correct answer. This test case is interesting, 
because a) the numerical results can be compared to analytical formulas and 
b) it is particularly difficult, in the sense that the dwell time, td ~ L/2vp, is 
(untypically) small (fp: Fermi velocity). To illustrate the method further, we 
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also apply it to a wire with four channels and show, that not only the con- 
ductance but also more complicated quantities, like the local current density, 
can be obtained. 



3.1 Single Channel Tight Binding Wire 



x-x-x+x-x+x-x+x-x 

L/2-2 V L/2 V L/2+2 V 

Fig. 5. Single channel tight binding wire with triple barrier realized by weak links 
(indicated by vertical double lines) as used in model calculations. 



Models and Analytical Results 

The model Hamiltonian of the clean tight binding wire is given by 

L 

H=-t/2Y,c\ci+r+c\^^Ci (12) 

i=l 

for spinless, non-interacting electrons. The corresponding dispersion relation 
reads 

e/s = — tcos(fca), (13) 
where a denotes the lattice constant, and for the density of states one has 

q{E) = , ^ (14) 

It exhibits the usual van Hove singularities at the edges of the band, which 
can also be seen in Fig. 6. 

In order to calculate the transport coefficient of the L— chain defined in 
(12) we should couple it adiabatically to a left and a right hand side reservoir. 
This can be done by attaching further half-infinite tight binding chains to the 
right and the left of the L— chain. The combined system is a perfect Id crystal. 
Its Bloch waves travel without any backscattering through the L— chain and 
therefore it is a perfect conductor with transmission unity: 

T{E)=l \E\ < t (15) 

Next, let us insert an obstacle, c. g. a strong triple barrier (sec Fig. 5), 
into the wire, a situation that still can be understood in all detail. The corre- 
sponding Hamiltonian H^^ is realized with hopping amplitudes occurring in 
(12) that take the values i=0.05 at the pairs of sites (i±=L/2,i±=L/2 ± 2) 
and t=l everywhere else. 
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Fig. 6. Clean single channel tight binding wire, upper panel: transmission. Ana- 
lytical result (15) (solid,thin) and numerical result for L=256 with fully absorbing 
boundary conditions (dashed; i(=32, /3=0.3, ?7=1). Also shown is numerical result 
with partially absorbing boundary conditions for L=16 (solid, wiggly line, is = l, 
P=oo, 77=1). Lower panel: density of states corresponding to the transmission curves 
shown in upper panel. 

The triple barrier has two eigenstates, a symmetric and an anti-symmetric 
one, which are energetically nearly degenerate since the center barrier is high. 
The energy is given approximately by the width of the double well inside the 
outer barriers, 3a. It corresponds to a wavenumber k — n/Sa which in turn 
implies a resonance energy close to t/2. Therefore, the transmission character- 
istics of the triple barrier should exhibit a superposition of two Lorentzians, 
one slightly below and one slightly above E = t/2. These are the features that 
can indeed be seen in the analytical result for the conductance (valid in the 
limit of weak coupling, V <S^t) 

where E±= cos(7r/3±y sin(7r/3)/3), t=l and 7=^^ sin(7r/3)3/24. A derivation 
can be found in appendix A. Here, we have reproduced for better clarity only 
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an expansion of the exact answer, given in (32), which is vaUd in the vicinity 
of the resonances. We also display the exact transmission in Fig. 7, upper 
panel. 




Fig. 7. Transmission of a single channel wire with a strong triple barrier (sec text). 
Upper panel: exact analytical result. Center: deviation of numerical results from 
analytical calculation: In t^^^^ (L=256, 1^=32, /3=0.3, r?=l, dashed; L=64, ij=16, 
(3=0.3, rj=l, dotted) with absorbing boundary conditions. Lower panel: partially 
absorbing boundary conditions (dotted, wiggly line: L=16, %=1, f3=oo, 77=1) (solid, 
smooth line: L=64, %=1, /3=oo, r)=l) 



Green's Function Method with Absorbing Boundary Conditions 

The transmission of the combined system - wire plus triple barrier - has been 
given already in (1). In the present case it reads 

g{E) = Tr,^ Gr.G^r^, (17) 

where the definition of G, (4), implies: 

G~^ = E — — — Sj^. (18) 
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The Hamiltonian H^^ of the extended molecule is given with Eq. (12) and 
the trace Tre^^ is over the corresponding Hilbert space. The operators r^ .jf^ 
are related to those pieces, Si,^^, of the self energy, S^^, which describe the 
level broadening due to the coupling to the left {l) and right {iC) leads: 

r,,^=i{s,,^-sl.^). (19) 

The precise form that the S^,^ take, depends on how we couple the wire to 
the external leads. In the spirit of section Sec. 2.3, we simply define as 
follows: 

Sij;_c = ivAj-, r]i=r]/{l + exp(3{i-is)) (20) 
Eij.^ = iriidij; r]i=r]/{l+ex.pl3{L-i-is)). (21) 

Three parameters have been introduced: 77 is the atomic Iciakage rate for all 
those atoms that are fully coupled to the outside; % describes the number of 
surface atoms on either side of the molecule; /3 is the adiabaticity parameter, 
that models a smooth transition into the external wire. 

Fig. 6, upper panel, displays the result of this procedure for the transmis- 
sion of the clean wire. As expected, the exact result (15) is recovered in the 
case with perfectly absorbing boundary conditions (abc) . For comparison, we 
also show a trace corresponding to incomplete absorption. The cavity modes 
manifest themselves in the transmission characteristics as relatively sharp res- 
onances. In order to highlight this aspect, the lower panel of Fig. 6 also shows 
the density of states of the wire. 

In Fig. 7 we present the transmission of the single channel wire with a 
triple barrier implantation. As can be seen from the upper panel, the agree- 
ment between the wire with fully abc and the analytical result is perfect. Once 
more we also display traces that result from a calculation with imperfectly ab- 
sorbing boundaries. Traces for two different cavity sizes, L=16, 64, are given. 
Like in the previous case, Fig. 6, the cavity eigenmodes give rise to system size 
dependencies, which are the spurious resonances in the transmission charac- 
teristics. By contrast, no remnant of the system size is left if perfect abc are 
used, see upper panel traces for _L=64, 256. 

Let us emphasize, that the good quality of the results Fig. 6 and 7 is not a 
consequence of fine-tuning parameters. We have ascertained, that the traces 
corresponding to perfectly abc are stable against variations at least in the 
parameter range i:=64-512, r?=0.1-10, /3=0.03-0.5, i^=8-64. 

In the test cases presented in this subsection, analytical results have been 
available in order to demonstrate that the choice of parameters associated 
with the absorbing boundary conditions was appropriate. In more realistic 
situations analytical results arc almost never available. Therefore additional 
criteria have to be given so as to establish that a certain choice of bound- 
ary conditions indeed provides sufficient absorption. The basic rule is, that 
a good implementation will yield transmission curves that are (largely) inde- 
pendent of the size N^^ of the cavity, of the choice of surface atoms inside 
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the cavity and of the atomic leakage rate rj. A calculation, that satisfies these 
requirements, is (usually) quite reliable. 

3.2 Local Currents in a Many Channel Tight Binding Wire 

As a further application of our method, wc calculate the local current density, 
jn, in a multi-channel wire. We begin by deriving a general formula relating 
to the Green's functions and self energies, that have been calculated in the 
preceeding section. Thereafter, we shall illustrate the result by calculating the 
local current distribution within a double well embedded in a four channel 
wire. 

Lattice Current Density in Terms of Green's Functions 
To start with, we consider the general model Hamiltonian 

H =^Y.*.y clc,'. (22) 

The multi-index u comprises the longitudinal, and transverse, £,y, wire 
coordinates. An expression for the local (longitudinal) current density may be 
obtained from the time dependent local density: 

= '^[HjC'^C/j] = — ^ i/ni/cjjCy — i/ii/ctc/n- (23) 

The component of the local particle current in the longitudinal direction 
(right), jp,, is given by the difference of those hopping events that enter site 
{ip,l^ from the left and leave there again: 

•^^~~2h invC^Ci, — t^j^clcp. (24) 

The expectation value of this operator is readily expressed in terms of the 
Green's function [13], G<^{t,t')=i{cl{t')c„{t)), 

In order to simplify the notation, we have introduced a name for sums like 
the one appearing in (24), which are restricted to the left/right half space: 
Since we only consider non-interacting fermions, G^^ takes a particu- 
larly simple form (e. g. [29]) 



G< =*G(/,r,+/^r^)Gt, 



(26) 
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which in the Umit of small voltages V leads to the following formula for the 
local charge current distribution = Jij-' 

e{3{i,))=-^V jdE ^ [Gir.-r.)G'U- (27) 

(■II V 

(When writing this expression, we have assumed for simplicity that H is time 
reversal invariant, so that t is real and t, G and F are symmetric matrices. 
Also, fc^x denote the Fermi-Dirac distribution of quasi-particles in the left 
and right hand side reservoirs.) 

Application: Double Dot Inside a Four Channel Wire 




-2-10 1 2 

Fig. 8. A double dot inside a tight binding wire consisting of four coupled strands 
(indicated by different symbols). 

In order to give an example for the usefulness of (27) , we calculate for 
a four channel wire with a double well barrier. The example illustrates, that 
the average current in the wire is a sum of contributions. They can undergo 
strong spatial fluctuations which are of the order of the mean current and in 
particular they can be positive or negative (backflow). It is only the sum of 
all of them which is independent of the longitudinal spatial coordinate. 

The clean 4-wire considered here is made up from four strands, which are 
the 1-wires given in (12). These 1-wires are arranged in a 4-fold cylindrical 
geometry and in transverse direction only nearest neighbors are being coupled, 
see Fig. 8. The coupling between all nearest neighbor pairs has been chosen 
t=2, except for those 9 pairs, that form the double well. They have t=0. 
The position of these "defects" can be given in terms of the longitudinal site 
index, i, and the transverse index ^, that labels the constituting 1-wircs in a 
clockwise fashion: ^ = 0, 1, 2, 3. We have switched off three couplings at and 
near the center {ic=L/2) of the wire, ictic ± 2, with site indices £=1,2,3. 
After the definition of the model Hamiltonian, JJeaf ) we also give the left and 
right contributions to the self energy , 
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Fig. 9. Local current density inside a four channel wire (I/=256, 1^=32, (5=0.1, ri=l) 
with a double dot (Fig. 8) at e=0.7 very close to a resonance of the transmission 
(see lower inset). Upper inset: behavior near the dots (position ±2, 0). Channels 1, 3 
(□, A, dashed line; same current density) and 2 (o,long dashed) are blocked, channel 
zero (o, dotted line) remains open. Inside the wells, j2 overshoots the transmis- 
sion (solid line) by 300 %, so that local backflow in the other current channels is 
generated. Main: sum of local currents (solid line) is conserved, i.e. independent of 
position (except for the surface region, where by adding the self energy leakage has 
been introduced). 

Eff.^ = 6ee'Sir,^, ^ = l,h (28) 

with a self energy per strand, Sij-x, as in (20, 21). 

Again, we may employ (18, 19), so that the numerical evaluation of (27) is 
straightforward (T^O). Fig. 9 shows the induced current density per applied 
voltage near the resonance energy, e = 0.7. Before we discuss this result in 
more detail, we first address the structure of the transmission function which 
is also displayed in Fig. 9, lower inset. It consists of three pairs of resonances, 
each pair resembling the symmetric and anti-symmetric eigenstate of the (iso- 
lated) double well. The pair of peaks closest to the band edges originates from 
hybridization of these states with those wire modes, that have a wavenumber 
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matching 7r/2a and an s-wave type symmetry in transverse direction. In these 
two peaks the current density is homogeneously distributed among the four 
constituting wires. The last sentence is not true for the remaining two pairs 
of resonances with p-wave character, where the current flows mainly in one of 
the wire pairs, cither (0,2) or (1,3). Evidently, the two resonances closest to 
e=l belong to the second category, since these are much sharper and less well 
split than all the others. 

Now, we can come back to the strong oscillations seen in the local current 
density of the resonance closest to the band center, main panel of Fig. 9. 
There, the current flow is mostly in the (0,2) pairs. Since due to the barrier 
these two channels are not symmetry equivalent, the current densities jio and 
ji2 can pick up different dependencies on the longitudinal wire coordinate, i. 
In fact, this must be the case since at the barrier position jo=l while j2=0. 

The phase locking between the local currents flowing along the (0,2) pair 
of strands has an interesting effect on the current flow inside each well: in this 
region, the component j2 acquires a value three times exceeding the average 
current flow. This value is compensated by a backflow in the other channels, 
so that a current vortex develops. 

4 Test Cases from Quantum Chemistry Calculations 

The purpose served by the tight binding calculations of the previous chapter 
was to demonstrate the principle. High precision in the calculations performed 
there was relatively easy to achieve, because the parameter space representing 
a perfect separation of energy (or time) scales was well accessible to numerical 
methods. 

In the case of quantum chemistry calculations that are feasible at present, 
the accessible system sizes often are not large enough in order to achieve such 
a clear scale separation. What we demonstrate in this section is, that our 
method operates reasonably well, also in a practical situation, where scale 
separation is not perfect. 

To this end, we shall consider two extreme C9/S6S, short gold chain with 
a conductance g ^ I and a molecule, di-thiophenyl with g <^ 1. Both objects 
are coupled to the tip of two tetragonal bipyramids of 14 Au-atoms each, that 
represent the extension modelling a piece of the electrodes, see Fig. 4. We 
define an effective single particle Hamiltonian, Hg^ , for these systems in the 
way explained in Sec. 2.1. 

The proposed construction mechanism has been formulated in coordinate 
(real) space. Therefore it matches particularly well with quantum chemistry 
calculations, that are using the local atomic basis sets \X,£) = \b) also intro- 
duced in Sec. 2.1. Adopting (5) to the present case, we introduce the following 
self energy: 

I]x = ivf^\b)[S-\b' {b'\, X = L,!!i. (29) 

5,6' 
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The overlap matrix, Shb' = {b\b'), that appears in this expression takes care 
of the fact, that basis states belonging to different atomic sites will not be 
orthogonal in general. ,^ is local, i.e. diagonal in the atomic basis set \b), 
in full analogy to (20, 21). Again the important input is in the strength and 
spatial modulation of the leakage function. In the present case, we choose 
it to be a constant, rj, for a subset of "surface atoms" and zero for all 
the others. In our calculations we take these sets to be the two layers of the 
pyramid (3x3 and 2x2) that are farthest from the molecule, see Fig. 4. 




E[eV] 

Fig. 10. Upper panel: Transmission of a two atom Au-chain linked to the tip 
of two tetragonal pyramids of 14 Au-atoms each. Self energy (29) has been used. 
Three traces with different rj are shown: 2.72 eV (dashed), 0.81 eV (solid), 0.27 eV 
(dotted). The D-symbol indicate the result with the canonical coupling and much 
larger pyramids with 55 atoms [3]. Lower panel: density of states calculated with 
the Green's functions used in the upper plot. Also the density of states for bulk gold 
is shown (symbol x). 
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4.1 Transmission of Au-Chain 

We begin our analysis with the two atom Au-chain. The transmission of such 
chains has been studied by various groups before [12, 30]. It is well known that 
reproducing the correct transmission curve is a sensitive test on the quality 
of the damping of the model. In Fig. 10, upper panel, we depict the 
transmission as obtained from the model (29). For comparison, also plotted 
is a result obtained from a much larger system with a self energy calculated 
directly from (3) (where has been replaced by the bulk Green's function 
Gs) [3]. Agreement for the two larger values of the coupling amongst each 
other and also with the original curve is established reasonably well with 
deviations typically less than 10 %. 

As was to be expected, the window of ry-values in which one finds good 
quantitative agreement between the various curves is relatively small. This is 
simply because of the very small cavity size. In fact, the "molecule" defined by 
the two Au-atoms in line, see Fig. 4, are separated from the effective surface 
regions 5l,r by only one gold atom. 

In the surface regions, the model self energy somewhat modifies the local 
material parameters, like the density of states (DoS). This can be seen in Fig. 
10, lower panel. The total DoS is strongly dominated by surface atoms. It has 
a dependency on 77 that is much stronger than the one of the transmission Fig. 
10, upper panel. Note, that this modification will have a substantial impact, 
if one were to set up a self consistent calculation with a local density obtained 
from the dressed Green's function (4). The present setup is not (and in fact 
does not need to be) self consistent in this sense, and therefore the modification 
of the surface spectral function will be without consequences for the transport 
calculations proposed here. 

We comment on the absence of a adiabaticity parameter /3 in the definition 
of the self energy (29) . The pyramids simulating the electrodes act as resonat- 
ing cavities in the same way that the single channel tight binding wire does, 
c.f. Sec. 3.1. However, the tight binding wire was special, in the sense that the 
number of surface atoms coupling to the infinite tight binding chain was only 
two, independent of the volume of the wire, L. For this reason the resonator 
modes had to be eliminated by introducing the adiabaticity parameter f3. In 
higher dimensions, the ratio of contact surface to volume is much more favor- 
able. Therefore the surface damping of the resonator modes is much stronger 
- no real need to introduce a /3-parameter here. 

4.2 Transmission of Di-Thiophenyl 

Finally, we apply our construction for the self energy to the paradigm of com- 
putational molecular electronics, the di-thiophenyl system. Again, the elec- 
trodes are modelled by the same pair of 14-Au pyramids, that we have used 
before for the 2- Au-chain. Accordingly, the construction of the model Hamil- 
tonian and, in particular, the self energy are just as in the previous section. 
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We will investigate two slightly different situations, where the sulphur 
atom, that tics the benzene to the Au-contact surface, once connects to a 
single Au-atom and once to three of them. 

S — Aui Coupling 




Fig. 11. Upper panel: transmission of di-thiophenyl with a S — Aui coupling (see 
Fig. 4). for different damping parameters rj: 0.27 eV (dotted), 2.7 eV (solid) and 27.2 
eV (dashed). Also shown 30- Au pyramids with conventional coupling (29). Lower 
panel: probability density P(e) for eigenvalues (proportional to the DoS). Three 
traces correspond to the transmission lines of upper panel. Symbols x indicate P{e) 
of bulk gold. 

The atomistic setup, that we consider in this subsection, is presented in 
Fig. 4. The S-atom acts as the barrier that disconnects the conjugated tt- 
system of the phenyl ring from the Au-atom forming the tip of the pyramid. 
This atom provides the separation for the junction from the contact region 
which is necessary in order to find results for the transmission independent of 
the choice of rj within a large parameter window. 

Indeed, our expectations are well confirmed by the numerical data. Fig. 
11, upper channel, shows transmission lines for rj varying over two orders of 
magnitude. All traces faithfully reproduce the salient resonance structures in 
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the vicinity of 1 eV about the Fermi energy, Ep = —5.05 eV. In this regime, 

the transmission is almost unaffected by the change in rj even though the 
average DoS changes by a factor of three, see Fig. 11, lower panel. Eventually, 
there is an impact of r] in the tails of the resonances, where the transmission 
is small and the DoS changes with 77 by an order of magnitude. This regime 
can be better controlled by creating an additional spacer layer of Au-atoms 
between the contact atoms and the junction. 

We believe, that the construction principle (29) is well suited to model 
the line broadening of resonances. There is, however, no prediction as for the 
line shift, of course. Line shifts occur, e. g. when charge reorganization takes 
place, involving charge flow from one subsystem to another. Such processes 
can be reliably modelled only by increasing the number of electrode atoms in 
the calculation and monitoring the results. 

S — AU3 Coupling 

In this paragraph, we investigate a situation in which the S-atoms couple to 
three Au-atoms rather than just one. There are two motivations for doing 
so. First, a two or three site configuration of the sulphur is energetically more 
favorable than the single site coordination used in Sec. 4.2 and thus more likely 
to be relevant for the understanding of experiments [3] . Second, we would like 
to present a simple example for a situation, where the dwell time in the cavity 
is not sufficiently long so that the parameter window has closed, in which the 
transmission traces are independent of the damping i]. 

In general, one avoids coupling the S-atoms directly to the surface layers, 
^L,R, because the change in the local DoS near the contact may feed back into 
the transmission. Therefore, for the S — Aus-coupling (Fig. 4) we include all 
Au-atoms into Sl.r except for those three Au-atoms, that bind to S. After 
this minor modification, calculations proceed in the same way as before. 

The resulting transmission is displayed in Fig. 12. For the 14- Au pyramids 
traces for three different damping values t; are shown. They vary appreciably 
one from another and the correct result is recovered only in relatively rough 
terms. The limited quality of the model coupling in the present case was to 
be expected. Since the number of surface atoms that couples to the S — Aua 
unit (eight) is much bigger than it was with the S — Aui coupling (four), also 
the dwell time is drastically reduced, too much for our simple model (5) to 
work with high precision. 

The situation is much more favorable if larger cavities are used. Fig. 12 
also shows the result for a 55-Au cluster which is in good agreement with the 
conventional calculation. 

5 Discussion and Outlook 

The method relying on absorbing boundary c;onditions (abc), that we have 
outlined in the preceeding sections, has an important advantage over the more 
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Fig. 12. Transmission of di-thiophonyl with a S — AU3 coupling (see Fig. 4). Thin 
traces: 14-Au pyramids with r/ 0.81 cV (dotted), 2.76 cV (solid) and 8.1 cV (dashed). 
Fat traces: 55-Au pyramids with r] 2.7 eV (dot-dashed) and conventional coupling 
(29) (solid). 

orthodox way of calculating an essentially exact self energy (3): it is much 
easier to handle. If the abc-calculation meets the requirements listed in Sec. 
(2.1), then both methods are expected to yield identical results for the Green's 
functions. 

We have demonstrated how the abc-approach can be combined with stan- 
dard quantum chemistry calculations for extended molecules, esW . Because it 
is not necessary any more to calculate self energies, , transport calcula- 
tions based on the Landauer-Biittiker theory are greatly simplified. 

The computational effort is dominated by the quantum chemistry calcu- 
lation for eiW . Since the number of electrode atoms, A^£ , that needs to be 
included in the calculation, is similar for both methods, one expects that also 
the computational effort is roughly the same. 

The number A^^ is too large in order to allow highly correlated, essentially 
exact methods to be used for calculating GeM ■ However, calculations based 
on density functional theory can be done very efficiently for these system 
sizes and also Hartree-Foek calculations are within reach. The last point is of 
particular interest, because this is a prerequisite to test functionals against 
each other (BP89, B3Lyp, LHF, etc.), that have a very different degree of self 
interactions. 

These tests are important for the recent debate on the origin of the dis- 
crepancy between theoretical and experimental results on the conductance 

observed for several organic molecules [2, 3, 4, 5, 6]. The conductance of 
monatomic chains can be investigated very well with the combination of meth- 
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ods, density functional theory (DFT) and Landauer-Biittiker formula, that 
have been employed in Sec. 4 [31]. By contrast, theoretical expectation for 
the transmission of organic molecules tend to deviate by one or more orders 
of magnitude from the experimental findings [3]. 

Apart from experimental difficulties also approximations, that are implicit 
to the DFT based transport formalism, could very well be responsible for this 
discrepancy. This is, roughly speaking, because the Green's functions, that 
derive from the (ground state) Kohn-Sham-formalism do not necessarily pro- 
vide a good description of the system dynamics. This description can indeed 
be acceptable, if electron density is relatively smooth and the system is at 
least close to metallic, like it is in single atom metal chains. Under these con- 
ditions, single particle wavefunctions, which are essentially plane wave states, 
give a good representation of the spatial properties of the true Green's func- 
tion. Then, it is mostly the (local) spectral properties (i.e. the band structure 
of the atomic wire) that determine the transport characteristics and those can 
be given quite accurately by present days DFT calculations. 

However, organic molecules arc a different case, because at least in the 
vicinity of the contact region, the electron density is far from homogenous. 
In these regions, the spatial structure of the KS-Green's function will in gen- 
eral not be a faithful representation of reality, certainly not within local or 
semi-local approximations of the exchange correlation functional. Since these 
very regions form exactly the bottleneck for the transport current, theoreti- 
cal calculations employing such functionals cannot necessarily be expected to 
be very precise. One can hope, that the useage of non-local functionals will 
improve upon this situation. Which of the various non-local terms, that exist 
in the exact (quasi-static non-equilibrium) functional, is the most important 
one and how to implement it in practical transport calculations, these are at 
present two of the most thrilling issues in the field of molecular electronics. 
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A A Triple Barrier in a Tight Binding Wire 

A one dimensional tight binding chain with only nearest neighbor hopping 
and three barriers is considered. 

x-x-x-x+x-x-x+x-x-x+x-x-x-x 

M-2 M-1 M Vl M+1 ... N Vj, N+1 ... O 0+1 0+2 0+3 

The Hamiltonian of this problem is: 

Hi J = ti + <5i+ij) (30) 

with hopping amplitudes ti given by 

tM = Vl left barrier 

In = Vc center barrier 

to = right barrier 

ti = 1 everywhere else 

and barriers located at positions M, N, O. The scattering states in the four 
different sections of the wire can be written as 

= e'''^ + re-'''^ left lead 

= ae**^^ + 6e~''^-' left to center barrier 

]Pj = ce*'^-' + de"*'^^ center to right barrier 

= te^^o right lead 

with energy E{k) = — 2cosfc (length scales are measured with respect to a). 
Let us now write down the consequences of the Schrodinger equation left and 
right of each barrier; 

E^M = ^M-l + Vl'^m+i 

EH/m+I = Vl^m + 
E^N = ^N~l + 

E^N+i = Vc'I'n + I'n+2 

E^o = ^o-i + VfltZ/o+i 
E^fQ+i = Vr^o + ^0+2 



or in matrix notation: 
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g(M+l)ifc g-(M+l)ife 

g(Af+l)ifc g-(Ar+l)ifc 
g(0+l)ife g-(0+l)ifeN^ 



y^g(M+l)ife y^g-(M+l)ife\ 



-Mik 



Vce 



e e 

(Ar+l)i/c -(]V+l)ifc 



te 



-,Nik 



(0+l)ik 



Vce- 

e 

Vr 
E - e^^ 



-Nik 



Since, our interest is in the transmission, T = we only extract an equation 
for the transmission coefficient, t: 



_(y^2g-(2M+l)«fe _ ^-{2M+l)ik-^^y^2^{2N+l)tk _ ^{2N+l)tkj'] (y^2gifc _ ^-tk^ 



2Ak 



2 -{2N+l)ik 



-{2N+l)ik 



.(y).2g-(2M+l)ifc _ ^-c^M+ljik^^Y^z^. 



-VLFcVRsin^fc. 



-ik 



2^{20+l)ik _ g(20+l)ife 



) 

(31) 



After specializing to the case of a strong, symmetric barrier, where Vl=Vc=Vr 
=V <C 1 and M=0, N^2, 0=4, we conclude 



T{k) = 2V- 



sin''(fca) 



I sin''(3fea) - sin^(2fca) + 2F2e*'^" sin(3fca) sin(2A;a) + o{V'^)\ 

(32) 

which is also displayed in Fig. 7, upper panel, in the body of the paper, where 
also a brief discussion of this result may be found. 
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